Draft version January 25, 2011 

Preprint typeset using K^'T^]X style emulatcapj v. 11/10/09 



ON THE ORTHO-.PARA RATIO OF H+ IN DIFFUSE MOLECULAR CLOUDS^ 
Kyle N. Crabtree-, Nick Indriolo'', Holger Kreckel^, Brian A. Tom^"', Benjamin J. McCall^"*'^ 

Draft version January 25, 2011 

ABSTRACT 

The excitation temperature Tgi derived from the relative intensities of the J — Q {para) and J — 1 
(ortho) rotational levels of H2 has been assumed to be an accurate measure of the kinetic temperature 
in interstellar environments. In diffuse molecular clouds, the average value of Tqi is ~70 K. However, 
the excitation temperature T{ll^) derived from the (J, K) = (1, 1) (para) and (1, 0) {ortho) rotational 
levels of H;^ has been observed to be '^30 K in the same types of environments. In this work, we 
present observations of Hjj" in three additional diffuse cloud sight lines for which H2 measurements 
are available, showing that in 4 of 5 cases Tqi and T(Hj^) are discrepant. We then examine the 
thermalization mechanisms for the ortho:para ratios of H^ and H2, concluding that indeed Tqi is 
an accurate measure of the cloud kinetic temperature, while the ortho:para ratio of H;^ need not 
be thermal. By constructing a steady-state chemical model taking into account the nuclear-spin- 
dependence of reactions involving HJ, we show that the ortho:para ratio of H;^ in diffuse molecular 
clouds is likely governed by a competition between dissociative recombination with electrons and 
thermalization via reactive collisions with H2 . 
Subject headings: astrochemistry - ISM 



1. INTRODUCTION 

Observations of H;^ in diffuse molecular clouds (diffuse 
clouds in which a significant fraction of the hydrogen is 
in molecular form: [Siiow &: McCallll200"6l ) have led to var- 
ious unexpected discoveries. The very first detection of 
H;^ along a diffuse molecular cloud sight line (toward Cyg 
0B2 12) showed an abundance a bout 10 times greate r 
than expected for the environment (|McCall et al.|[l998af ) . 
This surprising overabundance — also found toward sev- 
eral more diffuse cloud sight lines — led to the eventual 
conclusion that the ionization rate of H2 due to cosmic 
rays, C2, must be abo ut 1 order of magnitude larger than 
previ ously thought (|McCall et al.l 120031: llndriolo et"all 
l2007f l. Another puzzling outcome from H|^ observa- 
tions is that the average excitation temperature derived 
from the two lowest energy states ((T(H^)) « 30 K; 
llndriolo et al.l [2007') differs significantly from the aver- 
age excitation temperature derived f rom the two lowest 
rotational states of H p ((^0 1) ~ 70 K: lSavage et al.lll977l : 
iRachford et "all I200I [2009^ . Given that both species are 
expected to be thermalized to the cloud kinetic temper- 
ature by collisions, such a discrepancy is unexpected. 

Despite the fact that the J = and J = 1 levels of H2 
are essentially different "species" — conversion between 
the even-J para and odd-J ortho nuclear spin states 
is only possible through reactive collisions — it has long 
been assumed that Tqi is a good approximation for the 
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kinetic temperature, Tfei„, in diffuse molecular clouds. 
This is because H2 molecules are expected to experience 
many reactive collisions with protons during their life- 
times. As a result, the ortho and para populations of 
H2 should be brought into t hermal equilibrium w ith the 
proton kinetic temperature (|Dalgarno et al.iri973f) . 

Similarly, Hj^ also has ortho and para nuclear spin 
states {{J,K) = (1, 1) is the lowest lying para state, and 
{J,K) = (1,0) the lowest lying ortho state). For R'^ 
though, the relative population between the two states 
has been expected to be th ermalized by reactive coUi- 
sions with H9 (.McCall et al, 1998b; McCall et al . 200| 
IGibb et al.l[2O10l) . As the average values of Tqi and 
r(H^) do not agree, this does not seem to be the case. 

However, out of the 66 and 18 sight lines used to com- 
pute (Toi) and (T(HJ)), respectively, only 2 are shared 
between both samples. While previously reported val- 
ues of Tqi and T{'R^) do differ for t hese sight lines to- 
ward C Per and X Pe r bv about 30 K (ISavage et al.1ll977l : 
IRachford etHI 120021 : llndriolo et all |2007D, such a small 
sample does not provide particularly meaningful results. 
In order to increase the number of sight lines with or- 
tho and para column densities determined for both H2 
and , we have made observations searching for H;^ 
absorption features along sight lines with measured H2 
column densities. These observations and our data re- 
duction procedures are described in Section 2. Section 
3 discusses and justifies many of the assumptions made 
concerning molecular hydrogen and Tqi in diffuse molec- 
ular clouds. In Section 4, we examine the H;^ -I- H2 reac- 
tion and present steady state models to explore the tem- 
perature discrepancy between Tqi and T(Hj). Section 5 
reviews the observations and modeling, and presents our 
conclusions from the work. 



2. OBSERVATIONS & DATA REDUCTION 
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Target sight lines were se l ected based on H2 column 
densities (|Savage et al.lll977l:[Racliford et al.ll2002ll2009D 
and L-band magnitudes. We required that the J = 
and J = 1 column densities of H2 {N{0) and iV(l), re- 
spectively) both be known, and that the total H2 column 
density (iV(H2)) be greater than 10^° cm~^. The L-band 
magnitude was required to be brighter than 6 mag so 
that the necessary integration times would be relatively 
short. Observations focused on the 1)", i?(l,0), 
and i?(l,l)' transitions (at 3.668083 /xm, 3.668516 /im, 
and 3.715479 fim, respectively) which arise from the 
{J,K) = (1,1) and (1,0) levels of the ground vibra- 
tional state of H;^ , the only levels expected to be signif- 
icantly populated at average diffuse cloud temperatures 
(T ~ 70 K). 

Spectra in support of this pro ject were obtained 
using the Phoenix spectrometer (jHinkle et al.l I2003D 
at the Gemini South Telescope and the Cryo- 
genic High -resolution Infrar ed Echelle Spectrograph 
(CRIRES) (IKaufl et al.l [200 1) at the Very Large Tele- 
scope (VLT). Observations at Gemini South were made 
in queue mode, and the Phoenix spectrometer was used 
with its echelle grating and 0.17" slit to produce a re- 
solving power of about 70,000, and with the L2734 filter 
to select the order containing the it!(l, 1)" and i?(l,0) 
transitions. Observations at the VLT were made in ser- 
vice mode, and CRIRES was used with its 0.2" slit to 
provide a resolving power of about 100,000, and a refer- 
ence wavelength of 3715.0 nm to cover the R{1, 1)" and 
i?(l,0) transitions on detector 1 and the R{1, 1)' transi- 
tion on detector 3. The adaptive optics system was used 
with CRIRES to maximize starlight passing through the 
narrow slit. In addition to the science targets, bright, 
early-type stars were observed for use as telluric stan- 
dards. For all observations, the star was nodded along 
the slit in an ABBA pattern in order to facilitate the re- 
moval of atmospheric emission lines and dark current via 
the subtraction of neighboring images. A log containing 
the observed sight lines and respective integration times 
is given in Table [TJ 

Various standard IRAF^ procedures were used in the 
data reduction process. Given the different state of 
data available from Phoenix versus CRIRES, different 
amounts of processing were required for data from the 
2 telescopes. For each night of Phoenix observations, a 
bad pixel map was created from the average of several 
dark frames, and these pixels were interpolated over in 
the object and fiat frames. Flats were then combined 
to create a normalized flat field which was divided out 
of the object frames. Neighboring AB image pairs were 
subtracted from each other to remove atmospheric emis- 
sion and dark current. Finally, one-dimensional spectra 
were extracted using apall. In the case of CRIRES obser- 
vations, data were processed through the CRIRES spe- 
cific pipeline, resulting in fully reduced two-dimensional 
spectral images for each target (given the product codes 
SCOM and PCOM by the pipeline). One-dimensional 
spectra were extracted from these images. All such spec- 
tra from both telescopes were then imported to IGOR 
PrcQ where the remainder of reduction was performed. 

® http://iraf.noao.edu/ 
http://www.wavemetrics.com/ 
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Fig. 1. — Spectra showing absorption Hnes from the {J, K) = 
(1, 1) and (1, 0) states of Hg". Spectra for f Per and X Per were 
taken at U KIRT and are reproce ssed versions of the data previously 
reported in llndrioio et al] 120071) . The spectrum for HD 154368 was 
taken at Gemini South, while those for HD 73882 and HD 110432 
were taken at the VLT. Large deviations from flat continuum lev- 
els in the spectra for HD 154368, HD 73882, and HD 110432 are 
the combination of 2 effects: (1) artifacts due to continuum fitting 
at wavelengths far away from the H'^ absorption lines; (2) inabil- 
ity to remove the strong atmospheric methane line immediately 
shortward of the i?(l, 1)" line. Vertical lines mark the expected 
positions of the H^ absorption lines given previously determined 
interstellar gas velocities along each sight line. 

Individual Phoenix spectra within an exposure se- 
quence for a given target were then added together. Sci- 
ence target spectra from both telescopes were divided by 
telluric standard spectra to remove atmospheric absorp- 
tion features and to normalize each spectrum. These 
normalized spectra were wavelength calibrated with a 
typical accuracy of ~ 2 km s~^ using the vacuum wave- 
lengths of the atmospheric absorption lines. Calibrated 
spectra were then shifted into the local standard of rest 
(LSR) frame, and spectra of the same target from dif- 
ferent nights were combined using a variance-weighted 
mean. For each sight line the continuum surrounding 
the absorption features was then fit with an n*^ order 
polynomial (n = 7 — 9) and the spectrum was divided by 
the fit to re-normalize the continuum level. The result- 
ing spectra for sight lines with absorption features — 
HD 154368, HD 73882, and HD 110432— are shown in 
Figure [TJ Although H;^ spectra for ( Per and X Per 
(the other 2 sight lines wi t h bot h H2 and data) are 
reported in llndriolo et al.l (|2007l ). the reduction process 
utilized in that study did not combine spectra using a 
variance- weighted mean, nor did it fit fiuctuations in the 
continuum level with a polynomial function. For the pur- 
pose of consistency, we have reprocessed the data from 
both sight lines. No new data have been added, but 
differences in the method of processing have resulted in 
output spectra with slightly better signal-to-noise ratios 
(S/N). These spectra are also shown in Figure [TJ 

Absorption features due to Hj^ were fit with Gaussian 
functions in order to determine equivalent widths, veloc- 
ity full- width at half-maxima (FWHM), and interstellar 
gas velocities. Our fitting procedure uses the functional 
form of a Gaussian where area (as opposed to amplitude) 
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TABLE 1 

Observations 



Integration Time 



Object 


Date(s) of Observation 


Telescope 


(min) 


Standard 


HD 149404 


2009 Aug 05 


Gemini South 


32 


A Sco 


X Oph 


2009 Aug 30 


Gemini South 


12 


A Sco 


HD 152236 


2009 Aug 30 


Gemini South 


16 


A Sco 


HD 154368 


2009 Aug 30 


Gemini South 


32 


A Sco 


HD 53367 


2009 Dec 03 


VLT 


30 


K Ori 


HD 73882 


2009 Dec 03 


VLT 


30 


C Pup 




2010 Jan 17 


VLT 


30 


C Pup 


HD 110432 


2010 Jan 27 


VLT 


10 


Tj Cen 




2010 Feb 28 


VLT 


20 


r) Cen 




2010 Mar 02 


VLT 


40 


r) Cen 


fi Nor 


2010 Apr 05 


Gemini South 


72 


A Sco 



is a free parameter, and includes a fit to the continuum 
level, uq. In the case of the and i?(l,0) lines, 

both absorption features are fit simultaneously and a sin- 
gle best-fit continuum level is found. Uncertainties on 
the equivalent widths (SWx) and continuum level (Sy) — 
both at the la level — were output by the fitting process. 
To estimate the systematic uncertainties due to contin- 
uum placement, we forced the continuum level to yo + 6y 
and yo — Sy and re-fit the absorption lines. Variations 
in the equivalent widths due to this shift are small com- 
pared to those reported by the fitting procedure and so 
are not included in our analysis (i.e., a{W\)=6W\). As- 
suming optically thin absorption lines and t aking tran- 
sition dipole moments and wavelengths from iGoto et"al] 
(|2002f ) and references therein, column densities were de- 
rived from equivalent widths using the standard relation. 
All of these results are shown in Table [5] 

These observations increase the total number of sight 
lines with both and H2 detections from 2 to 5. Col- 
umn densities, para-fractions, and excitation tempera- 
tures for both species along all 5 si ght lines are co l lected 
in Table |3l H2 data com e from I Savage et all ()1977f ) 
and IRachford et all (j2002D . Uncertainties on all values 
are la. The excitation temperatures inferred from the 
i?(l,0), and absorption lines of H+ 

range from 20 < r(H;^)< 46 K, while those reported 
for H2 vary from 51 < T(H2) < 68 K. In 4 sight lines Tqi 
is greater than r(H;^) by about 30 K, while for X Per 
Toi and T(H+) are consistent within uncertainties. Still, 
these observations clearly show that for the same inter- 
stellar material along 4 different diffuse molecular cloud 
sight lines the excitation temperatures derived from Hj^ 
and II2 do not agree. 

3. H2 THERMALIZATION 

Given the large discrepancies between Toi inferred 
from II2 and the excitation temperature of , it is im- 
portant to re-examine the assumption that the inferred 
Toi accurately reflects the kinetic temperature of the dif- 
fuse molecular clouds. If this assumption were incorrect, 
it would be conceivable that H;^ provides the true (lower) 
temperature of diffuse molecular clouds. There are at 
least 4 reasons this assumption could be invalid: (1) ob- 
servational errors in the determination of the J = and 
J = 1 column densities of H2; (2) an insufficient fre- 
quency of H+ -|- H2 collisions to achieve steady state; (3) 
the steady state of this reaction being different from the 
thermodynamic equilibrium; and (4) errors caused by a 



varying J = : 1 ratio along the line of sight. In the 
following subsections, we investigate each of these possi- 
bilities in turn. 

3.1. Observational Determination of H2 Columns 

The measurement of the column densities of J = and 
J = 1 of II2 is performed by profile fitting to spectra of 
the Lyman (A-X) band in the ultraviolet, recorded with 
Copernicus or FUSE. The II2 transitions are optically 
thick, and are completely opaque in the line cores. Given 
the difficulties in accurately retrieving column densities 
from optically thick transitions, one might worry that 
the inferred Tqi is contaminated by uncertainties caused 
by radiative transfer. 

According to B. L. Rachford (private communication, 
2010), the detailed shape of the combined profile of the 
J = line and the two J — 1 lines is quite sensitive to 
the ratio of the column densities of these two rotational 
levels, and thus provides a very useful probe of Tqi. Be- 
cause multiple vibronic bands of H2, which are known 
to be relatively free of contamination by stellar lines, 
are used in the determination of Tqi, it is difficult to 
envision any systematic errors that could affect the mea- 
surements. The statistical errors in the column density 
measurements are typically ~0.1 dex, and we can there- 
fore conclude that the ultraviolet measurements provide 
an accurate and fairly precise measure of the ratio of 
A^(0) to iV(l) for H2. 

3.2. Frequency of + H2 Reactions 

The J = and J — 1 rotational levels of II2 belong 
to different spin modifications {para and ortho, respec- 
tively) and are therefore not thermalized by non-reactive 
collisions or radiative transitions. Only chemical reac- 
tions, in which the protons are exchanged, can affect the 
nuclear spin modification of II2EI 

Since H2 is formed in a highly exothermic reaction 
on grain surfaces, its nascent rotational (and spin) dis- 
tribution is presu med to represent a high temperature 
(iTakahashil 120011) . In the high temperature limit, II2 
should be formed with an ortho:para ratio of 3:1. If 
an insufficient number of reactive collisions occur be- 
tween the formation of an H2 molecule and its subsequent 

^ Strictly speaking, large inhomogeneous magnetic fields, such as 
found on the surface of paramagnetic catalysts, can also affect the 
nuclear spin modification of H2, but we assume that such effects 
are negligible in interstellar conditions. 
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TABLE 2 
Absorption Line Parameters 







fLSR 


FWHM 






N{J, K) 


a{N) 


Object 


Transition 


(km s^-"-) 


(km s^-"-) 


(lO-f* ^m) 


(10-6 tim) 


(10" cm-2) 


(1013 cm-2) 


Q Per 




7.7 


11.0 


0.99 


0.13 


4.09 


0.53 




R{1,0) 
R{l,l)" 


6.1 


9.0 


1.00 


0.11 


2.53 


0.29 


X Per 


8.2 


9.1 


0.80 


0.17 


3.34 


0.69 




R{1,0) 
_R(1,1)" 


6.3 


10.2 


1.30 


0.18 


3.29 


0.45 


HD 154368 


5.4 


6.0 


1.79 


0.30 


7.43 


1.24 




R{1,0) 
R{l,ir 


5.2 


5.6 


1.12 


0.29 


2.83 


0.74 


HD 73882 


5.9 


3.9 


1.44 


0.21 


5.97 


0.86 




R{1,0) 


5.7 


3.2 


1.16 


0.19 


2.94 


0.48 




i?(l, 1)" 


5.4 


3.5 


1.34 


0.15 


6.15 


0.69 


HD 110432 


-3.8 


6.9 


0.74 


0.06 


3.08 


0.24 




i?(l,0) 


-3.3 


7.5 


0.83 


0.07 


2.11 


0.17 




i?(l,l)' 


-3.1 


8.1 


0.69 


0.06 


3.15 


0.28 



Note. — Column 3 (ijlsr) gives the interstellar gas velocity in the local standard of rest frame. Column 
4 (FWHM) gives the full width at half-maximum of the absorption features. Columns 5 and 6 show the 
equivalent width, W\, and its la uncertainty, (t(Wa), respectively. Columns 7 and 8 give the column 
density of H^ in the state each transition probes, N{J,K), and its uncertainty, cr{N), respectively. Values 
for these parameters in the ( Per and X Per sight lines were previously reported in Indriolo ct al. ( 2007). 
The new values for both absorption lines toward Per and the R{1, 0) line toward X Per are consistent with 
the previously published results within uncertainties. However, the new and old results for the -R(l, 1)" 
line toward X Per are inconsistent. Upon inspection, we found this to be due to a bad fit to that line 
during the 2007 ana lysis. In all cases, the values published herein should be taken to supersede those from 
llndriolo et all l(2007l 1. 



TABLE 3 
Target Sight Line Properties 







C Per'''^ 


X Per^'^ 


HD 154368== 


HD 73882== 


HD 110432'='^ 


HJ Results 


Ar(l,l) 
iV(l,0) 


(10^3 cm-2) 
(10^3 cm-2) 

(K) 


4.09 ±0.53 
2.53 ±0.29 
0.62 ±0.04 
28 ±4 


3.34 ±0.69 
3.29 ± 0.45 
0.50 ±0.06 

46i?J 


7.43 ± 1.24 
2.83 ±0.74 
0.72 ± 0.06 
20 ±4 


6.08 ±0.12 
2.94 ±0.48 
0.67 ±0.04 
23 ±3 


3.11 ±0.05 
2.11 ±0.17 
0.60 ±0.02 
30 ± 2 



H2 Results 



log[Ar(0)] 
log[Af(l)] 



(cm 
(cm 



(K) 



20.51 ±0.09 
20.18 ±0.09 
0.68 ±0.06 
58 ±6 



20.76 ± 0.03 
20.42 ±0.06 
0.69 ±0.04 
57 ±4 



21.04 ±0.05 
20.54 ±0.15 
0.76 ±0.07 

51 ±8 



20.99 ±0.08 
20.50 ±0.07 
0.76 ±0.05 
51 ±6 



20.40 ± 0.03 
20.27 ±0.04 
0.57 ±0.03 

68 ±5 



Note. — Measured column densities for the lowest lying ortho and para states of H2 and Hg are shown 
for the 5 sight lines with all such data available. Also shown are the para-fractions for each species and 
the rotational temperatures derived from a simple 2-state system analysis. 

^ Updated analysis of Hg" data originally presented in llndriolo et all 120071 ) 
^ H2 data fromLSavage et al.l 019771) 
" H2 data from Rachford ct al. ( 2002) 

May be affected by multiple velocity components ()Crawfordlll995l ) 
« P-H;^ fraction: N(l, 1)/{N{1, 0) + N{1, 1)) 
f P-H2 fraction: N{1)/{N{0) + N{1)) 



destruction (by photodissociation or cosmic-ray ioniza- 
tion), then the average ortho:para ratio of H2 may he 
somewhere between the nascent value (3:1) and the ther- 
mahzed value (1:2 at 60 K, for example). This could lead 
to Tqi overestimating the true cloud kinetic temperature. 

The number of reactive collisions suffered by an av- 
erage H2 molecule in its lifetime can be expressed as 
^frxn = Tiife/Trxn, whcrc Tuff, is the average hfetime 
of an H2 molecule and Trxn is the average time between 
reactive collisions. If A/'rxn 3> 1, then the ortho:para ra- 



tio of H2 should reflect the steady-state of the reaction 
in question. 

In diffuse molecular clouds, II2 is formed on grains 
at a rate of Rn}in{}l), where R is the grain formation 
rate (typical l y taken to be ab out 3 x 10^^'' cm"^ s~^; e.g. 
ISpitzedllQTStiGrv et al.l[200l . nn = n(H)-|-2n(H2) is the 
total number density of hydrogen nuclei, and n(II) is the 
number density of atomic hydrogen. The destruction of 
II2 is dominated by cosmic-ray ionization and photodis- 
sociation (following absorption in the Lyman bands) , and 
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has a rate of (C2 +r)n(H2), where F is the photodissocia- 
tion rate. In steady state, these two rates are equal, and 
we can solve for Tufe = (C2 + F)"^ = n(R2)/[RnHn{l{)]. 
By using the definition of the local molecular fraction 
/^^ = 2n(H2)/riH (which we will simply call /), we can 
rewrite this expression as rn/e = //[2i?nH(l — /)]• 

Assuming that the reaction of H+ + H2 dominates the 
interconversion of 0-H2 and p-H2(or</io- and para-H2, 
respectively), we can write Trxn = [fcic'T-(H+)]~^, where 
kic is the rate coefficient for the interconversion reaction. 
Substituting into the equation for Mrxm we find 

_ fc.e n(H+) / 
R nn 2(1-/)- 

To estimate the number density of H''", we consider 
the steady state of its formation and destruction. H+ 
is formed by cosmic-ray ionization of H atoms, at a 
rate of Ch"-(H), where i s the cosmic-ray ionizatio n 
rate of H (2.3Ch ~ 1-^2] 'Glassg old fc Langeri HqtI . 
Given the abundances of various species in diffuse clouds 
and the rate coefficients for reacti ons between H+ and 
such species, chemical models (e.g. IWoodall et al.|[2007[ ) 
predict that H+ is destroyed predominantly by charge 
transfer to atomic oxygen, with a rate of kctn{R^)n{0). 
Equating these rates of formation and destruction and 
solving for n(H+), we find n(H+) = Ch(1 — f) /[kctx{0)], 
where a;(0) = ri(0)/nH. Finally, substituting this into 
the expression for Mrxn gives 

^ kic Ch / 1 
" R kct 2nu x{0) ■ 

Adop ting va l ues of kic = 2.2 x 10~^° cm'^ 
s-i (IGerlichI [TOOOh . R = 3 x IQ-^^ cm^ s-\ 
C2 = 4 X 10-1^ s-i (llndriolo et al.l 120071) . 
kct = 7.31 X 10-i°(r/300)°-^3g-22 5.9/T ^^3 

(IWoodall et al.l[200l. x(q ) « 3 x lO"'' (|Cartledge et al.l 
120041 iJensen et al.l l2005| ). and typic al diffuse cloud 
value s of / = 0.^| and nn = 100 cm'^ (|Snow fc McCalll 
[2006I) . we find that at T - 70 K, Kxri 1400. Thus, the 
typical H2 molecule will experience over 1000 reactive 
collisions during its lifetime, and we can safely assume 
that the initial orthoipara ratio of H2 is irrelevant; 
instead, the observed ortho:para ratio should reflect the 
steady state of the reactive collisions. 

3.3. Steady State of + H2 Reactions 

From a thermodynamic perspective, one would ex- 
pect that the steady state of the H+ -I- H2 reaction 
would represent a thermal distribution of 0-H2 and p- 
H2, if no other processes influence the spin modifl- 
cations. This expectation has been confirmed by a 
phase space theoretical calculation by IGerlichI ()1990[ ). 
who found that the ortho:para ratio could be ex- 
pressed at low temperatures (30-80 K) by the expression 
9.35 exp(— 169. 4/T), quite close to the thermodynamic 
expectation of 9 exp(— 170.4/r). Evidently the nuclear 

^ While the Hne-of-sight molecular fraction is lower, that quan- 
tity integrates over atomic gas not associated with the diffuse 
molecular cloud. We are therefore using the local molecular frac- 
tion typical of diffuse molecular clouds. 
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Fig. 2. — The line-of-sight rotational temperature, Tqi, is de- 
termined for a cloud containing both a warm and cold component. 
The cold, inner component is varied between 10 K and 100 K, while 
the warm, outer component is set at 100 K. Different curves show 
the results for cases where the warm component contains 1/2, 1/4, 
1/8, and 1/16 of the total material in the cloud, and are labeled 
accordingly. Cross hairs mark the inferred inner cloud values of 
Tqi given the observed line-of-sight values of Tqi and estimated 
fractions of H2 in a 100 K cloud component for the 5 sight lines 
considered herein. From left to right the cross hairs mark HD 
73882, HD 154368, C Per, X Per, and HD 110432. 

spin selection rules for this chemical reaction, which sup- 
press the ortho-to-para conversion by a factor of 6, do 
not significantly impact the final distribution. 

3.4. Line of Sight Integration Effects 

One remaining concern regards the estimation of Tqi 
in a diffuse molecular cloud from the column densities of 
J = and J = 1, which are integrated quantities along 
the line of sight. If some of the II2 resides in hotter, 
mostly atomic gas where H;^ is not abundant, that hot 
H2 would cause the observed line-of-sight Tqi to exceed 
Tqi in the molecular cloud. We expect that such con- 
tamination would not be a major effect, as II2 is known 
to self-shield very effectively from the interstellar radi- 
ation field; thus, the amount of II2 in primarily atomic 
(and presumably warmer) gas is likely to be quite small 
compared to the amount of H2 in the molecular cloud 
itself. 

To estimate the magnitude of this effect more quan- 
titatively, we take a simple cloud model with a hotter 
outer region and cooler inner region. Assuming that 
Tqi — 100 K in the outer region (bas ed on Copernicus ob - 
servations of diffuse atomic clouds: Ijenkins et al.lll983l ). 
we varied Tqi in the inner region between 10 K and 100 K 
for a set of models where the outer region contained 1/2, 
1/4, 1/8, and 1/16 of the material in the cloud. We then 
computed the line-of-sight Tqi that would be derived con- 
sidering both regions of gas. The result of this analysis 
is shown in Figure [21 

Taking the diffuse cloud model of INeufeld et al.l (|2005D 
with a constant density (nn = 100 cm~^) and standard 
UV radiation field (xuv = 1) illuminating the cloud 
from one side, we then determined the amount of II2 
expected to be in the outer region for comparison with 
observed II2 column densities in diffuse clouds. We de- 
fine the transition from the outer to inner regions to be 
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at E{B -V) = 0.04 (A^h « 2.3 x lO^o cm-^), about 
half of the total color excess (hydrogen col umn density) 
found to supply self -shieldinR effects in H2 (jSavage et al.1 
[19771 iGillmon fc Sh ull 2006). Integrating n(H2) in the 
outer region gives Nouter{^2) = 6 x 10^^ cm~^. Because 
this model effectively only treats one side of a cloud, we 
compare this value to half of the total H2 column densi- 
ties reported in Table [H The two extreme cases are HD 
110432 and HD 154368, where the outer region accounts 
for about 1/4 and 1/12 of the total cloud, respectively. 
Taking the observed values of Tqi and using the appro- 
priate curves on Figure 2, we estimate the temperature 
of the inner cloud region for each of the 5 sight lines 
considered herein. The results are marked in Figure 2 
as cross hairs, and show that the line-of-sight values of 
Tqi overestimate the inferred inner cloud values of Tqi 
by only about 5-10 K. As such, the observed values of 
Tqi should be relatively close to the true values of Tqi in 
molecular clouds. We therefore assume for the remainder 
of this paper that the line of sight Tqi does represent the 
diffuse molecular cloud kinetic temperature. 

3.5. Summary on H2 Temperature 

From the preceding discussions, we conclude that UV 
measurements accurately measure the column densities 
of J = and J = 1 of H2 in diffuse clouds, that the chem- 
ical reaction between H+ and H2 occurs ~ 10'^ times 
during the life of an average H2, and that the steady 
state of this chemical reaction leads to an ortho:para ra- 
tio that closely reflects the kinetic temperature of the 
gas. Furthermore, we conclude that it is unlikely that 
the integration along the line of sight introduces signif- 
icant contamination of the inferred Tqi. In summary, 
then, measured values of Tqi in diffuse molecular clouds 
should accurately reflect the cloud kinetic temperature. 
Consequently, the excitation temperature of , which 
is signiflcantly lower than Tqi, must not always reflect 
the kinetic temperature. 

4. ORTHO AND PARA H+ 

Having shown that the temperature discrepancy be- 
tween Tqi and T(K^) in diffuse molecular clouds is real 
and that Tqi accurately reflects the cloud kinetic temper- 
ature, we now consider the processes related to H^ ther- 
malization in diffuse molecular clouds to examine why 
T(H;^) might not match the kinetic temperature. 

4.1. Nascent and Thermalized Para-H^ Fractions 

Hi^ is produced by the reaction 

H2+H+->H++H, (1) 

which follows the cosmic-ray ionization of H2 to form 
HJ. The nuclear spin modiflcation of the product H;^ 
depends on the nuclear spin modiflcations of the reactant 
H2 aii d H^ according to the selection rules given bv lOkal 
(I2OOI . It is most convenient to express the nuclear spin 
modifications in terms of the para-fractions, rather than 
the ortho:para ratios, so we define 

^ n(p-H2) 

" n{p-R2) + nio-H^) 

and 

n{p-R+) + n{o-R+y 



As the cosmic-ray ionization of H2 is not expected to 
affect the nuclear spin modification, we can further as- 
sume that the para-fraction of H J is also given by p2 ■ Ta- 
blelUdemonstrates, using these reactant fractions and the 
nuclear spin branching fractions, that the para-fraction 
of newly formed H3 is ps = (1/3) + (2/3)p2, assuming 
that the rate for the Hj" -I- H2 reaction is independent of 
nuclear spin configuration. 

In diffuse molecular clouds, the vast majority of the 
H2 population lies in the lowest ortho and para states, 
as the temperature of 70 K is well below the energy of 
the next states (the J = 2 state lies 510 K above J — 0, 
and J = 3 lies 844 K above J — 1). Therefore, to de- 
rive p2 from astronomical observations we use the for- 
mula p2 = A^(0)/[iV(0) + N{1)]. H+ on the other hand 
does have energetically accessible para states, as the (2,2) 
and (2,1) states lie only 151.3 and 249.2 K above the 
(1,1) ground state. However, the (2,2) and (2,1) states 
are expected to quickly undergo radiative decay to the 
(1,1) state at the temperatures and densities of the dif- 
fuse molecular clouds (Oka & Epp 2004). Furthermore, 
population in the next ortho st ate, (3,3) has n ot been 
observed in these environments (iOka et al.ll2005() . so the 
vast majority of o-H^ is in the (1,0) state. Consequently, 
to calculate p^ from the astronomical observations, we 
use p3 = iV(l, l)/[iV(l, 1) + iV(l, 0)]. 

Figure [3] shows the nascent pa distribution as a func- 
tion of p2. This figure also shows the total para-fraction 
of a thermalized sample of H;^ at various tempera- 
tures, calculated using t he energy levels E{J, K) from 
iLindsav fc McCalll (|200lD . Also plotted in the fi gure are 
the results of the astronomical observations presented in 
Table[3l In diffuse molecular clouds, pa generally appears 
to lie between the nascent pa and the thermal pa values, 
suggesting an incomplete thermalization of the nuclear 
spin modifications o f H j" . 

As discussed by lOka fc E"pdI (|2004D . the aforemen- 
tioned spontaneous emission from the (2,2) and (2,1) 
states decreases T(H^) relative to Tqi. They show that 
for a cloud density of 100 cm-^ and 60 < Tqi < 120 
K, T(H+) should fafl in the range of 40-50 K, and this 
accounts for a ~ 40 — 80 K difference between Tqi and 
T(H;^). However, in the temperature ranges discussed 
here, the pora-fractions of H^ and H2 are nonlinear with 
respect to these excitation temperatures. In terms of 
Pa, all T(Hj^) above about 40 K should have about the 
same pa, while pa changes substantially when T(Hjj') falls 
below 40 K, as can be seen in Figure [3l Spontaneous 
emission will raise the apparent pa [as derived from the 
N(1,1):N(1,0) ratio] relative to the thermalized pa and 
consequently lower r(Hj) with respect to Tqi. However, 
Figure [3] illustrates that the magnitude of this effect can- 
not account for the discrepancy observed in the astro- 
nomical observations in these environments in terms of 
the para-fractions, with the possible exception of X Per. 

4.2. The Reaction of with H2 

As in the case of H2 , the nuclear spin modifications of 
H;^ cannot effectively be changed by radiative transitions 
or by non-reactive collisions; only chemical reactions can 
do so. In this case, the reaction H;^ -I- H2 — >■ (H^)* H2 
-|- H;^ is the most efficient mechanism for interconverting 
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TABLE 4 

Nascent Para Fraction 



Reaction Collision Fraction Branching Fraction p-Hif Fraction 





h P-H2 


{P2)2 


1 


Pi 




p-hJ - 

o-H^ - 
o-H^ - 


h 0-H2 


P2(l - P2) 


2/3 


(2/3)(l - 


P2)P2 


h P-H2 


(1 - P2)p2 


2/3 


{2/3)(l - 


P2)P2 


h 0-H2 


(1-P2)2 


1/3 


(1/3){1 - 


-P2? 



Total 



(1/3) + (2/3)p2 



Note. — This table presents the calculation of the nascent p-H^ fraction 
formed in diffuse molecular clouds from the + H2 reaction, assuming 
that cosmic ray ionization of H2 to form does not affect its nuclear spin 
configuration. The collision fraction represents the fraction of total + 
If2 collisions with the specified nuclear spin configurations. The branching 
fractions are for p- Hj" f ormation, and are derived from nuclear spin selection 
rules ( Quacl3 |l977l : [OE a 2004). The final column presents the calculation of 
the nascent p-H^ fraction. 
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Fig. 3. — The nascent (dashed line) and thermalized (numbers, 
in K) p-H!^ fraction as a function of the P-H2 fraction. The circles 
represent the observations of diffuse molecular clouds summarized 
in Table [3l with Icr uncertainties. 

o-H^ and p-H^ ■ When H;^ and H2 collide, there are 
three possible reaction outcomes: 



H3+ 
H+ 



H 



2^ HJ+H2 (identity), 
H2 ^ H2H+ + H2 (hop), and 
H2 ^ HH+ + HH (exchange) 



(2) 
(3) 
(4) 



In the case of reaction [5J the nuclear spin configura- 
tions of the and H2 remain unchanged, while in re- 
actions [3] and S] the configuration may change. However, 
like reaction [1] the hop and exchange pathways must 
obey nuclear spin selection rules (Quack 1977; Oka 2004). 
For instance, in order for a reaction between p-H'^ and p- 
H2 to form o-Hj^ the reaction must be an exchange, and 
0-H2 must also be formed to conserve the total nuclear 
spin angular momentum. 

A potential energy surface based on high-le vel ab ini- 
tio ca lculations is available for the system (|Xie et al.l 
120051) . Based on the surface stationary points, a bar- 
rier of 52.2 cm~^ must be overcome for a hop reaction 
([3]) to occur, and a barrier of 1565.9 cm~^ for an ex- 
change reaction (|4]) to occur. The dissociation energy 



is calculated to be 2903 cm ^; therefore (H^)* formed 
from association of with H2 has sufficient energy to 
overcome these barriers. As such, the complex may un- 
dergo many hop and exchange processes over its lifetime, 
and given sufficient time, the product distribution may 
approach a statistical limit. As the statistical weights 
for the hop and exchange reactions are 3 and 6, respec- 
tively, the branching ratio a = f^t^ is 0.5 in the sta- 
tistical limit. Quantum reactive scattering calculations 
are presently unavailable on this potential energy sur- 
face, so experimental data are necessary for determining 
the value of a, as well as its temperature dependence. 

The only experimental deter mination of a for t he Hj " 
-I- H2 system was performed bv lCordonnier et all (j2000t) . 
This study was done by spectroscopically measuring the 
ortho-to-para ratio of Hj^ formed in a discharge of pure 
P-H2 at 400 K, and under these conditions, the value a — 
2.4 was derived. No measurements at lower temperatures 
have been reported for this system, but the isotopically- 
substituted reaction + H2 has been studied at a va- 
riety of collision energies using an ion trap/guided beam 
technique ( Gerlich 1993). It was observed that a varies 
substantially with the Dj^-H2 collision energy. As this en- 
ergy decreases, a approaches the statistical value of 0.5, 
and the value 2.4 is reached at an energy correspond- 
ing to the average collision ene rgy at ^440 K, in general 
agreement with the study by iCordonnier et all |2000) . 
However, a direct comparison of these results to Hj^ in 
the interstellar medium is problematic owing to the en- 
dothermicity of the reaction channel and the nonthermal 
reactant internal state distribution in the experimental 
measurement. 

The final consideration for this reaction is the frac- 
tion of reactions that lead to no change in the nuclear 
spin modification, S***^. A large value for S'*'^ would in- 
dicate that nuclear-spin-changing collisions are a small 
fraction of the total number of H;^ + H2 collisions, and 
the thermalization process would be slower than the colli- 
sion rate. In fact, there is experimental evidence for this, 
as studies of the H^ + HD H2D+ + H2 reaction give 
a rat e coefficient of 3.5 x 10"^" cm"^ s^^ (jGerlich et al.l 
l2002f ). much slower than the Langevin rate coefficient 
These results lead to S'-'^ ^ 0.8, but 
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it is possible that S'^'^ could be different for the purely 
hydrogenic system, which is tliermoneutral rather than 
exothermic. 

4.3. Steady State Para-H^ Fraction from + H2: 
"Bimolecular Reactive Equilibrium" 

After taking into account its chemical physics, does the 
steady state of the + H2 chemical reaction lead to 
a completely thermalized in the interstellar medium? 
To consider this question, we have constructed a simple 
steady-state model for ortho and para-H^ , in terms of 
nuclear-spin-dependent rate coefficients kxxxx for each 
potential sub-reaction (e.g., koppo- o-H^ -|- P-H2 — p- 
+ 0-H2). The derivation of this model, which we 
shall call the "bimolecular reactive equilibrium" (BRE) 
model, is presented in the Appendix, and the resulting 
expression for is given in the Appendix as equation 
(iBl) . 

The rate coeffic ients themselv e s wer e computed using 
the prescription of Park fc Lightl ()2007 t) , which takes into 
account both the nuclear spin branching fractions as well 
as energetic considerations at the state-to-state level, us- 
ing a microca nonical approach. This work has since been 
extended by iHugo et all (|2009[ ) to deuterated versions 
of this chemical reaction, and the latter authors report 
quantitative agreement between the two sets of calcula- 
tions. We therefore judge these rate coefficients to be 
reliable within the context of this theoretical approach. 

The Park & Light code (provided by K. Park, private 
communication 2009) requires five input parameters: the 
kinetic temperature, the rotational temperature of 
and H2, and the three branching fractions S^'^, S'^°p, and 
gexch^ For these calculations, the rotational temperature 
was held at 10 K in each nuclear spin manifold in order 
to have the vast majority of o-Hj]' in (1,0) and p-H'^ in 
(1,1). Therefore, we express the inputs to the model in 
terms of only three parameters: T = Tkin, S^'^, and a, 
as S"^ + S''°P + S""""^ = 1 and a = S^°p / S""""^ . The 
code then outputs all of the rate coefficients required in 
equation (jB4p . For a single set of branching fraction 
values, the rate coefficients were calculated for 10 < T < 
160 K in steps of 10 K, and p2 was set to its thermal 
value for each calculation. 

Figure |4] shows the results of the BRE model for a 
fixed value of a = 0.5 for various values of S^'^ ranging 
from 0.1-0.9; similarly. Figure [5] shows results for fixed 
gid _ Q ^ g^j^j various a ranging from O-oo. The re- 
sults of the calculation are not particularly sensitive to 
the fraction of collisions that are reactive (as traced by 
S'*'^) or the ratio of the hop to exchange outcomes (a). 
Since in all cases p^ falls near its thermal value, these 
results suggest that regardless of the values of a or 5*'', 
the H;^ -I- H2 reaction should essentially thermalize the 
nuclear spin species. This stands in clear contradic- 
tion to the reported astronomical observations in diffuse 
molecular clouds, with the exception of X Per. The dis- 
crepancy between Tqi and T(H^) cannot be explained by 
the BRE model, and must then be explained by a lack 
of equilibration via this chemical reaction. 

An interesting aspect of these results is that the steady 
state ps at some temperatures is actually below the value 
of 0.5 expected based on statistical weights alone (often 
called the "high-temperature limit"). This appears to 
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Fig. 4. — BRE calculations of the para-H'^ fraction as a function 
of the para-H2 fraction, under the influence of the Hg" -|- H2 reac- 
tion. The thin solid line represents the thermal limit (as in Figure 
(3)1, and the circled crosses represent the results of our calculations 
(based on Park &; Light's model) for a = 0.5 and various values 
of S^"^ ranging from 0.1 (small purple) to 0.9 (large brown). Each 
cluster of crosses represents a calculation at a single temperature, 
ranging from 10 K (upper right) to 160 K (lower left) in steps of 
10 K. Also plotted are the nascent Hg" distribution and the astro- 
nomical observations. 
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Fig. 5. — Same as Figure|4l except that 5*'' is held at 0.5 and a 
varies between (small purple) and 00 (large brown). 



be a robust result for a > 0.5, at least in the range of 
S^'^ explored here. This effect may have been observed 
experimentally in our group's recent measurem ents of 
supers onically expanding hydrogen plasmas. iTom et al.l 
(|2010f l reported ps = 0.491 ± 0.024 for a normal hydro- 
:en (p 2 = 0.25) plasma at ~ 80 K, and iKreckel et al.l 
20101) reported p3 = 0.479 ± 0.02 in a warmer (- 200 
K) normal hydrogen and argon p lasmaP^ finally, in 



work outside our laboratory, Kreck el et al.l ([2007) have 
reported P3 = 0.4 in a low-temperature ion trap. All of 
these measurements suggest that it is in fact possible to 
achieve p^ < 0.5, and lend some evidence to support the 
calculated results. 

^'^ Recent measurements in our laboratory have confirmed, with 
higher statistical significance, p3 < 0.5 in hollow cathode plasmas 
containing normal hydrogen. These results will be reported else- 
where, but it should be noted that p3 in these plasmas may be 
influenced by three-body collisions due to the higher number den- 
sities. 
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lision rate for + 
sions an will experience is then 



4.4. Steady State Para-H^ Fraction from + H2 and 

H+ +e- 

We now consider whether there are enough reactive 
collisions within the lifetime of an average H;^ in dif- 
fuse molecular clouds to bring the spin modifications 
into BRE. The destruction of in such clouds is domi- 
nated by dissociative recombination (DR) with electrons, 
and the lifetime is simply Vife = (fcz?flTi(e ))~^ (the 
reciprocal of the destruction rate), where kon is the 
dissociative recombination rate coefficient. The aver- 
age time between reactive collisions can be expressed 
as Trxn = (fcrc?i(H2)) ~^ , whcre krc is the reactive col- 
H2. The average number of colli- 

■ '^lifej ^rxn — 

[krc/kDR\ [n(H2)/n(e~)]. 

Assuming for the moment that koB. is independent of 
the nuclear spin modification, we adopt a typical value 
(for T - 70 K) of koR = 2 x 10"^ cm^ s"! (|McCall et al.l 
120041) . The ratio n(H2)/n(e~) can be rewritten as f /2xe, 
where Xe is the electron fraction, typically 1.5 x 10~^ 
assuming cha rge neutrality and that C"*" is the domi nant 
ionic species (ICardelh et all [1991 ISofia et al.l[200l . If 
we adopt / = 0.9, and take krc to be the full collision 
rate of H+ + H2 (1.5 x 10"^ cm^ s"!; lAdams fc Si^ 
Il987| ). we find that Nrxn ~ 20. However, if we instead 
adopt the s maller reactive rate coefficient ^ 3 x 10"^" 
cm^ g-i of iGerlich et al.l (|2002D . we find that Nrxn ^ 5. 
With such a small number of collisions in the lifetime of 
Hj~ , p3 may not reach the value predicted by equation 
IB4I In the appendix we de rive a more complete steady 
state expression (equation IC7p including the effects of 
both the H;^ -t- II2 reaction and nuclear-spin-dependent 

DR rates (fce,p and fce,o for p-^^ and o-\i^). 

We call this model simply the "steady state" model, 
and we adopt the values / = 0.9 and = 1.5 x 10~ 
as before. Figure [6] shows the results of this steady 
state model if we assume that the DR rate coeffi- 
cient is the same for both nuclear spin modifications 
(we have adopted the temperature-dependent value of 
iMcCall et ai1[200l . In this case, the values of depend 
quite sensitively on S^'^, as this represents the fraction of 
-I- H2 collisions that are nonreactive during the rel- 
atively short lifetime of H;^. Consequently with higher 
values of S*"* (larger circles in Figure [6]), the p-H;^ frac- 
tion in steady state is closer to the nascent fraction. For 
S'*'' = 0.9, which corresponds to a reactive rate coeffi- 
cient of krc = 1-9 X 10~^° cm^ s~^, the calculated are 
in reasonable agreement with most of the observations. 
The upper range of the X Per uncertainty is consistent 
with S"^ = 0.7. 

In Figure [7l we instead consider the calculated DR 
rate coefficients for ortho and para-\i^ presented in 
Idos Santos et al.l (|2007l ). Their prediction is that p-Hj]" 
is destroyed considerably faster by electrons at low tem- 
peratures than o-W^ ; consequently, even for large values 
of S'*'', the steady state pz does not approach either the 
nascent or astronomically-observed values. In the ab- 
sence of the Hj]' -I- H2 reaction {S^'^ — 1), pz would be 
governed by a steady state determined by the competi- 
tion between the formation and the spin-dependent DR 
processes, and this is shown in Figure |8l 
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Fig. 6. — Steady state calculations of the para-Hg fraction as a 
function of the para-H2 fraction, under the influence of both the 
+ H2 reaction and dissociative recombination. The plotted 
quantities are analogous to those in Figure |4] In this case, the o- 
Hg" and p-Hg" DR rate coefficients fce^o and fce^p are assumed to be 
equal and taken from lMcCall et al.l I I2004I ). Each vertical cluster of 
points represents a calculation at a single temperature, beginning 
at 160 K in the lower left and decreasing by 10 K each point moving 
to the right. 
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Fig. 7. — Same as Figure [B] except using th e spin-dependent 
dissoc iative recombination rate coefficients from Idos Santos et al.l 
pOOTi 'l. 
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Fig. 8. — Steady state calculations of the para-Yl^ fraction as 
a function of the para-Y{2 fraction, under the influence of forma- 
tion and dissociative recombination only. The solid line shows the 
thermal limit, and the dotted line represents the results of our 
calculations, where we have used the spin-dependent dissociative 
recombination rate coefficients from dos Santos ct al. (2007). Also 
plotted are the nascent Hj}" distribution (dashed line) and the as- 
tronomical observations. 
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If th e calculated rate coefBcients of Idos Santos et all 
(|2007f ) are correct, it is difficult to explain the observed 
Ps- This is because, with the exception of X Per, the 
observed are higher than the curve resulting from 
the steady state of H^^ formation and destruction us- 
ing these DR rate coefficients, and inclusion of the H;^ 
+ H2 reaction further drives ps toward the value ex- 
pected for th ermal equihbrium . Re cent storage r i ng ex - 
periments bv lTom et all (|2009[ ) and lKreckel et all (|2005[ ) 
both saw an increased DR cross-section when is pro- 
duced fro m iP-H2; however, re cent imaging results pre- 
sented in iKreckel et al.l ()2010D suggest that the ions in 
these experiments have been heated during extraction 
from the ion sources, and the difference between the o- 
H;^ and p-ii^ may therefore have been overestimated. 
Further experimental work is clearly needed to pin down 
the enhancement (if any) in p-H^ DR, and confirmation 
of the theoretical predictions would also be quite helpful. 

To summarize, according to our models the reaction of 
Hjj' with H2 is expected to effectively thermalize the nu- 
clear spin configurations of H;^ at steady-state, provided 
that sufficient collisions occur within the lifetime of an 
Hj^. In diffuse molecular clouds, however, the average 
number of reactive collisions with H2 suffered by an H;^ 
is small, indicating that the formation and destruction 
rates of the two nuclear spin species should be impor- 
tant. A more complete model which takes these factors 
into account reaches reasonable agreement with observa- 
tions in 4 of 5 sight lines provided S^'^ is on the order 
of 0.9 and o-H^ and p-H^ are destroyed at equal rates 
owing to DR. Reconciling the observations with the sp in- 
dependent theoretical rates of Idos Santos et all ()2007[ ) is 
difficult, and accurate experimental measurements of the 
spin-dependent DR rates of Hg at low temperature are 
needed. 

5. CONCLUSIONS 

While all evidence seems to suggest that Tqi inferred 
from ultraviolet spectroscopy of H2 accurately reflects 
the kinetic temperature of diffuse molecular clouds, the 
observed excitation temperature of II;|^ is clearly non- 
thermal in 4 of the 5 mea sured sight lin es. Based on the 
microcanonical model of iPark fc Lighti ([2007,), we have 
constructed a steady state model to predict the para-H^ 
fraction (pa) if reactive collisions between and H2 
control the spin modifications of Hj^. Those results show 
P3 slightly below the limit expected for full thermaliza- 
tion, and far from the observations. However, a steady 
state model that incorporates both the Hj^ -I- H2 reac- 
tion as well as the formation (following cosmic-ray 
ionization) and destruction (by electron recombination) 
can reproduce the observed ps if the reactive collision 
rate is somewhat slow and the dissociative recombination 



rates for ortho and para Hg are comparable. Our inter- 
pretation, given the currently available data, is that 
suffers relatively few spin-changing collisions with H2 in 
its lifetime, and is thus incompletely equilibrated by this 
reaction. The observed para-H^ fraction therefore lies 
between the nascent fraction and the nearly-thermal frac- 
tion that would be reached with sufficient reactive col- 
lisions. If our model is correct (and the spin-dependent 
DR rates of are nearly equal at low temperature), 
this marks the first determination of the reactive rate 
coefficient of the H;^ -I- H2 reaction, and suggests a value 
on the order of 10^^" cm^'^ s^^. 

Fully quantum reactive scattering calculations of the 
Hj^ + II2 reaction would be highly desirable, as they 
would pin down the state-to-state rate coefficients needed 
to predict the interstellar para-H^ fraction. Further ex- 
periments and theoretical calculations to elucidate the 
dependence (if any) of the dissociative recombination on 
the nuclear spin modification of H;^ are also urgently 
needed. Once the effects of the reactive collisions and dis- 
sociative recombination are fully understood, the para- 
fraction in diffuse molecular clouds can be calculated 
as a function of the kinetic temperature and the ratio 
of the molecular fraction to the electron fraction. This, 
in turn, suggests that Hj^ may become a useful "ther- 
mometer" for diffuse molecular clouds with high extinc- 
tion, where ultraviolet measurements of H2 are not fea- 
sible. However, the calibration of this thermometer will 
require further experimental and theoretical efforts. 
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APPENDIX 
REACTIONS AND RATES 

In this appendix, we derive the formulas for the bimolecular reactive equilibrium and steady-state para-H^ fractions. 
We consider 4 processes: cosmic ray ionization of H2, formation of , the Hj^ -I- H2 reaction, and dissociative 
recombination of Hj^. If all nuclear spin configurations are considered, this gives a total of 28 reactions. The chemical 
reactions used in the model and their rate expressions are summarized in Table [5l It should be noted that for this 
section, we employ the chemist's notation of using square brackets to refer to the number density of the respective 
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TABLE 5 

Reactions and Rate Equations Used in Models 



Number Reaction Rate Comments 



1 


Ha + 


CR ^ H+ + 

H+ + 


e- + CR' 


C2[Ha] 




Cosmic ray ionization 


2 


Ha + 


H 


fci[H2][H+J 




Hg" formation (see Table [TJ 


3 


i-Hj 


+ i-Ha m- 


■H+ + n-Ha 




'][i-Ha] 


Thermalization reaction for Hg" 


4 


0-H+ 


+ e" ^- Ha 


+ H (or) 3H 


fce,p[p-H+]| 
fce,c[0-H+]| 




p-Hf DR 
o-nj DR 


5 


+ e" ^ Ha 


+ H (or) 3H 





species. 



Note. — The branching fractions for o-Hj and p-Hj in reaction 2 are assumed to be 
given by nuclear spin statistics. In reaction 3, i, j, m, and n represent the nuclear spin 
configurations of the respective species (o or p). Some of these 16 reactions are forbidden by 
nuclear spin selection rules, and others are not used directly in the derivation because they 
do not result in a change in the Hg" nuclear spin configuration. Square brackets refer to the 
number density of the species. 



DERIVATION: BIMOLECULAR REACTIVE EQUILIBRIUM MODEL 



Consider the case that Hjj" formation and destruction are slow compared with the + H2 reaction. We can then 
ignore formation and destruction processes and write the rate equation for p-H^ only in terms of the latter reaction: 

^[p-H;^] = {{koopo + koopp) [0-H2] + {koppo + koppp) b-H2]} [o-H;^] 

) [0-H2] + {kppoo + kppop) [P-H2]} [p-H;^]. (Bl) 

Assuming steady state, equation |BT] is equal to 0. The right side can then be divided by [Hg ][H2] in order to express 
the rate in terms of p2 and as follows: 

= {{koopo + koopp) (1 - P2) + [koppo + koppp) P2} (1 - Pa) 

- {{kpooo + kpoop) (1 -P2) + {kppoo + kppop) P2}P3- (B2) 

The resultant equation can be solved for p^: 

^ {koopp ^" koopo){^ P2) ~1~ {koppo ~H koppp)p2 (B3) 

{koopp ~t~ koopo ^~ kpoop ^ kpooo) {^ P2) ^ {koppo ^" koppp ^" kppoo ^" kppop)p2 

Owing to nuclear spin selection rules, the rate coefficients koppp and kppop are rigorously 0. Removal of these terms 
gives the final result: 

^ {koopp + koopo){^ P2) "I" koppoP2 (B4) 

{koopp ^" koopo H~ kpoop ~\~ kpooo){^ P2) ~t~ {koppo ~t~ kppoo)P2 
DERIVATION: STEADY STATE MODEL 

Consider now the case in which H;^ formation and DR compete effectively with the H;^ + H2 thermalization reaction. 
We make the assumption that formation of p-^l^ from + H2 is governed only by the nuclear spin branching fractions 
presented in Tabled Under these conditions, the full rate equation for p-H^is: 

^b-H+] =fci(b-H2][p-H2+] + ^b-H2][o-H+] + ^[o-H2]b-H+] + i[o-H2][o-H+]) 

) [0-H2] + {koppo + koppp) [p-li2]} [0-H3 ] 
{(^pooo H~ kpoop) [o-i^2] H~ {kppoo H~ kppop 
~k,Je-][p-ll+]. (CI) 

From TablelU we can reduce the first line in this equation to fci[H2][H^]{(l/3) + (2/3)p2}- We now invoke steady 
state arguments for all species. For HJ, C2P2] = A:i[H2][HJ], therefore fci[H^] can be replaced by ^2- Setting the 
resultant equation equal to zero and dividing by [Hj]'] [H2] as before gives: 
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+ {{koopo + koopp) (1 - P2) + {koppo + koppp) P2} (1 - Ps) 
~ {{kpooo + kpoop) (1 ~ P2) + (kppoo + kppop) P2} pz 



(C2) 



This equation can be further simphfied by examining the steady state value of [H||], which begins with the equation 
C2P2] = fce[H|^][e~]. More specifically, if we include the possibility for different DR rates for o-H^ and p-Hj^, we 
obtain the equation: 



C2P2] = [e-](fce,p[p-H+] +fc,,„[o-H+]). 
Dividing both sides of equation IC3l bv [H3 ][H2] results in an expression for C2/[H3 ] 

C2 



-{ke,pP3 + fee, 0(1 -Ps))- 



Substituting this relation into equation IC2I gives 





= ||^(^e,pP3 + fce.o(l -P3))(^ + ^P2) 

+ {{koopo + koopp) (1 - P2) + (koppo + koppp) P2} (1 - Ps) 
{{kpooo ~t~ kpoop ) (1 -P2) + {kppoo + kppop)p2}p3 



TPs- 



Solving for and removing the koppp and kppop terms yields: 



(C3) 



(C4) 



(C5) 



P3 



ke,o Ijfj] i,\ ^" fP^) + {koopp + koopo){^ ~ P2) ~\~ koppoP2 



ke,p [1^2] (3 3^-^) ke,o \\^2] (3 ^ 3^-^) ^ {koopp ^ koopo ^" kpoop ^" ^pooo)(l P2) ^ {koppo ~t" kppoo)p2 

Finally, the ratio [e^]/[H2] can be replaced by 2xe//, which results in equation IC7I 



(C6) 



P3 



ke,o^y^ + §^2) + {koopp + fcoopo)(l ^ P2) + koppoP2 



ke.p^^ (I " 1^2) + fce,o% (I + IP2) + (A:o 



+ ^pooo)(l ~ P2) + {koppo + kppoo)P2 



(C7) 
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